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Abstract 

Puzzling properties of sodium cobaltates Na^CoO^ become an issue 
of numerous recent theoretical and experimental investigations because 
of very high thermopower and intricate phase diagram of this mate- 
rial. Experiments indicated Na 2 /3Co02 to have a kagome-like cobalt 
sublattice, which implies additional complexity of this strongly cor- 
related system. In the present work we have employed the ab-initio 
GGA+U modeling method to investigate spin states of Na2/3Co02, 
using an exact crystal cell established in the experiment. Stability of 
sodium imprint for concentration x = 2/3 proposed in the publica- 
tions has been obtained. An AFM-A type magnetic ordering has been 
recognized as a ground state. Moreover, a competition between low 
spin state (LS) and intermediate spin state (IS) has been established 
to depend on on-site Hubbard U correction parameter. LS state has 
been identified ground state of Na 2 /3Co02- 
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1 Introduction 



The intricate ground-state properties of the sodium cobaltates Na^CoC^ 
are a subject of many recent theoretical and experimental investigations. 
This material holds much promise for thermoelectronics due to its large 
thermopower [1] together with the relatively low resistivity [2]. Magnetic 
fields strongly affect the thermopower in this material, which suggests that 
electron-electron strong correlations and spin degrees of freedom are impor- 
tant. The discovery of superconductivity with T c about 5 K in Na .33CoO 2 • 
1.3 H 2 [3] revived the interest to the physics of the lamellar sodium cobal- 
tates providing the first known example of a non-cuprate oxide superconduc- 
tors. 

Such properties as the superconductivity [3], high thermopower coeffi- 
cient [2] , uncertainty of the charge state of the cobalt ions in triangular crystal 
lattice and strong influence of the sodium ions arrangement onto electronic 
properties of C0O2 planes [4] make the investigation of sodium cobaltates 
very important for the modern physics of condensed matter. 

Structurally Na^CoC^ consists of quasi two dimensional C0O2 layers 
stacked along the crystallographic c-axis, separated by Na layers. Form- 
ing the triangular lattice, the planar Co ions are surrounded by edge-sharing 
06 octahedrally, which are trigonally distorted along c. 

The global T-x phase diagram of sodium cobaltate Na^OoC^ [5] has two 
regimes - with ferromagnetic and antiferromagnetic correlations. In addition, 
there are narrow insulating state at x — 1/2, which separates two distinct 
metallic regions, as well as superconducting state at x 0.3. A phase tran- 
sition between two regimes of magnetic correlations appears at x* m 0.65 at 
low temperatures. At x < x* the spin dynamics display predominantly an- 
tiferromagnetic correlations at low T, contrary to the ferromagnetic ones at 
x > x*. An origin of this magnetic transition remains a challenging problem, 
which motivates the present study. 

From the experiments it is known that phase above x* ~ 0.65 corresponds 
to long range AFM order of A- type, i.e., spins are ordered ferromagnetically 
in C0O2 plane, whereas exchange integral between planes is antiferromag- 
netic. It is supposed, that spin and charge ordering of Co 3+ /Co 4+ at x — 2/3 
are due to the simultaneous Na + ordering and magnetic transition in C0O2 
plane in metallic state. But the mechanisms of a such complex transition 
still remains the open problem. 

Minor deviations of the local electrostatic potential on cobalt sites, such 
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as those induced by local Na + order or disorder, may therefore significantly 
influence the electronic state [4]. Thus inplane sodium ordering is in a great 
importance. In Ref. [6-8] sodium ordering patterns for different sodium con- 
centration (x = 0.5, 0.6, 2/3, 0.75) were proposed. To predict the ground 
state, authors have used Monte Carlo simulations with GGA cluster expan- 
sion method. But in bulk samples not only inplane ordering, but also the 
relative position of sodium planes is very important. In the next section 
influence of the sodium ordering will be discussed in more details. 

1.1 Electronic structure 

In the cobaltate the electronic 3d level of the Co ion is splitted by a cubic 
crystal field of oxygen octahedron into the lower t 2g triplet and higher e g 
doublet orbital state with a crystal gap about 2 eV [9]. A trigonal distortion 
of 6 octahedra further splits the t 2g triplet into the doublet e' g and singlet 
dig state. In a Co 3+ valence state with d 6 electronic configuration the six 
lower energy levels fog) are filled, with a total spin S = 0, while Co 4+ should 
only retain one hole in the t 2g multiplet, with S = 1/2. Experimentally, 
differentiation between the valencies of cobalt ions can be induced by the 
sodium ordering [10]. It was suggested by Alloul et al [10] that such spatial 
arrangement of Co 3+ /Co 4+ valence state is induced by a specific periodic 
distribution of Na + ions, because nearby sodium cation could attract the 
electron onto the cobalt and fulfill the lower t 2g triplet. However, the micro- 
scopic origin of such ordering of the Co valence state and its relation with 
the simultaneous magnetic transition was not still established. 

Hereafter we will consider different spin configurations of cobalt ions Co 34 " 
and Co 4+ , notations can be found in the Table 1. There ZS stands for zero 
spin, HS for high spin, whereas LS stands for low spin and IS for intermediate 
spin states. 

In the present work we investigate an appearance of the unique state at 
the triangular Co lattice in the cobaltates at x = 2/3: the metal with the 
robust ferromagnetic spin and charge (valence) arrangement. To this end 
we study a role of the coexisting realistic mechanisms: (i) the d-electronic 
interaction with a non-uniform electric potential induced by the superlat- 
tice of the Na + ions in Na 2 / 3 Co0 2 and (ii) the strong Coulomb repulsion 
of the electrons in the 3d shells. There are several methods to deal with it 
- Hubbard model, DMFT and other. We treat the electronic correlations 
within the Hubbard model and employ an ab-initio calculations within the 
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GGA+U (LSDA+U) approximation to obtain the states of the system. Ab 
initio calculations with Hubbard model for strong correlation (LSDA+U) of 
simplified two dimensional model lead to the absence of small pockets at the 
Fermi surface [11, 12]. But the simplified models used in these works - for 
example, a single C0O2 layer, doped with electrons - do not demonstrate 
some peculiarities of the electronic ground state, for example, a wide band 
gap for x = 1/2. Thus, using a realistic crystal cell as well as taking into 
account strong correlations of d-electrons is in a great importance. 

Table 1: Spin configurations for cobalt ions Co 3+ and Co 4+ 



Co 3 + 


ZS 




HS 


configuration 


l 2g 




+5 1 

l 2g e g 


spin 


S= 





S = l 


Co 4 + 


LS 




IS 


configuration 


l 2g 


e° 


+4 I 

l 2g e g 


spin 


S= 


1/2 


S=3/2 



2 Ab initio modeling methods 

First-principles calculations were performed by means of the spin-polarized 
DFT+U method with Hubbard correction. Core electron states were rep- 
resented by the projector augmented-wave method [13] as implemented in 
the Vienna Ab Initio Simulation Package (VASP 5.2.2) [14], which is em- 
bedded into MedeA® interface. The Perdew - Burke - Ernzerhof exchange 
correlation revised for solids (PBEsol, see Ref. [15]) and a plane- wave rep- 
resentation for the wave function with a cutoff of 500 eV were used. Both 
internal coordinates and unit-cell lattice parameters were fully relaxed. The 
Brillouin zones were sampled with a 5 x 5 x 5 mesh including the gamma 
point. The density of the mesh for all calculations is approximately one 
point per 0.004 A~ 3 . SCF convergence criteria was set to 10~ 5 eV. To deal 
with strong correlations the simplified (rotationally invariant) approach to 
the LSDA+U, introduced by Dudarev et al. [16] was used. The Hubbard 
U c s = U — J value in the Hamiltonian is taken to be in general 5 eV for Co 
(but for certain cases we provide an analysis for range of U). This value for 
U is between the values of U — 4.91 eV for Co 3+ and U = 5.37 eV for Co 4+ 
obtained with first-principles perturbation theory in LL r Co0 2 [17]. While 
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GGA+U (or LDA+U) is most often used to open up a Hubbard gap in the 
electronic structure, it was previously shown that in related systems it also 
has a significant and meaningful effect on the energy. As it removes the 
self-interaction on the d orbitals, it leads to a strong charge localization and 
significant changes in the ground-state structures [18]. 

During our calculations we have tested two PAW potential for sodium: 
with seven valence electrons (2p 6 3s 1 ) and with one valence electron (3s 1 ). It 
was found that, in the first case, all p-electrons are localized on the sodium 
cation and do not change spin and charge distribution for cobalt ions. Thus 
the usage of sodium potential with only one valence electron is sufficient. 

3 Calculation model 

By means of NMR/NQR techniques, which are very useful for dealing with 
local properties of ions, several unique crystallographic positions of sodium 
and cobalt ions was identified and minimal crystallographic cell with sym- 
metry group R — 3c was built based on this information (Figure 1) [10]. 



Figure 1: Crystallographic cell of Na2/3CoO"2, deduced from NMR/NQR ex- 
periments. See [10] for more details. 

As it can be seen from Figure 1 the two general crystallographic posi- 
tions of cobalt ions can be distinguished. First position Col (which has two 
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types - Cola and Colb) has the minimal distance to the sodium ion. Atom 
in Col positions form a 2D charge ordered state - a honeycomb network of 
Co 3+ ions with t\ g e° g electronic configuration. Cobalt ions in another crys- 
tallographic positions Co2 (Co2a and Co2b) are forming a kagome lattice of 
atoms with average 3.44+ valence state. It means that the electrons are not 
fully attached to its own ions, but they are hopping along kagome lattice. 
Calculation shows that in a\ g and e' g orbitals the cobalt electronic states are 
hybridized with the oxygen one. 

Schematic electronic configurations for one symmetry irreducible layer of 
Na 2 /3Co02 in low spin and intermediate spin states are depicted on Figure 2. 

(a) (b) 

3Col 3+ 5Co2 3+ 4Co2 4+ 3 Col 3+ 5 Co2 3+ 4 Co2 4+ 




a lg -H--t a lg -W— *— 

it it it it it 

ZS ZS LS ZS HS LS 

Figure 2: Electronic configuration diagram of one symmetry irreducible layer 
of Na2/3Co02 in low spin(a) and intermediate spin (b) states. The total 
magnetic moment of layer is 4 /ie and 14 /ie, correspondingly. 

We perform geometric optimization of elementary cell in both the LS 
and IS spin states within GGA and GGA+U methods, within range of U 
from 1 to 6 eV. It should be noticed, that even without symmetrization of 
charge densities and forces, the arrangement of sodium ions remains the same. 
Thus, the stability of proposed sodium spatial distribution is confirmed. The 
energies of system under consideration in FM and AFM interlayer order in 
LS state have been calculated. The energies comparison indicates that the 
AFM order has less energy of about 7.5 meV per one irreducible layer and 
thus is more favorable. 

The charge/spin separation for optimized structures in LS and IS states 
was compared. We start from the case of the absence of strong correlation 
between 3d electrons (the Hubbard parameter is equal to U = 0). Then 
the GGA provides a uniform with FM spin-polarized metallic LS state of 
Co-planes. In turn, increasing of Hubbard repulsion U parameter results in 
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the localization of magnetic moments at Co2 sublattice when U exceeds ps 3 
eV. The spin state is obtained to be non-uniform: in the LS state one finds 
the robust magnetic moment of about 0.6 /is at the Co2a(b) sites, forming 
the kagome superlattice, while Col sites remain almost spinless. Thus a 
kind of a "spin separation" between the sites in the triangular Co lattice is 
appeared. Such a type of separation resembles the spin separation indicated 
by the spin-polarized GGA+U calculations for Nao.5Co02 [19]. 

For the IS state the situation is more complicated. It is not possible to 
stabilize IS state for the low values of Hubbard U parameter. The minimal 
value of Hubbard parameter at which the IS state was obtained is U = 4 
eV. This value is close to the U — 5 eV reported in [20], which is necessary 
to stabilize high spin state in Nao.6iCo02- So, one can maintain that our 
calculation is coherent with the other ab-initio calculations. 

Spin separation curves for IS state are depicted on Figure 3,(b). One can 
conclude that the increasing of Hubbard U parameter leads to the growth 
of magnetic moment at cobalt ions in "favorable" crystallographic positions 
(Co2b and/or Co2a), i.e. to the localizations of holes in Co2 sublattice. 

For the geometric parameters, such as the lattice constants a and c, there 
is no single motif. Increasing of U leads to decreasing of hexagonal inplane 
lattice parameter a in LS state, but small growth in IS state (Figure 4, (a)), 
whereas in the direction which is perpendicular to Co0 2 planes the lattice 
constant c demonstrates monotonic growth for both spin states (Figure 4, 
(b)). The maximal relative difference of lattice constants between two spin 
states is 1.5% for a parameter and 0.9 % for c parameter, which lies in the 
usual DFT error range of 2-3 %. Probably, one can explain such behavior 
by redistribution of electron charge density from inplane toward interplane 
orientation, which leads to repulsion of layers. Significant increasing of both 
lattice parameters in IS state can be explained in the following way: in the 
IS state electrons on cobalt ions occupy the upper e g bands. These bands 
are geometrically oriented toward nearest oxygen ions, which has partial 
negative charge and, as results, tends to repulse from cobalt ions, which 
leads to increasing of both a and c lattice parameters. 

Experimental lattice parameters for sodium concentration x close to 2/3 
are presented in the Table 2. As it can be seen, experimental lattice param- 
eters are closer to IS state, but the relative difference of order 2 % lies in 
usual range of geometrical parameters error, calculated within DFT meth- 
ods. Thus, the lattice parameter can not be a strict criteria to identification 
of experimental spin state. 
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Figure 3: Magnetic moments of cobalt ions depends on on-site Hubbard U 
correction parameter in low (a) and intermediate (b) spin states. 
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Figure 4: Lattice parameters a (a) and c (b) in a dependence on on-site 
Hubbard U correction parameter in LS and IS states. "GGA" denotes the 
pure GGA method, without LSDA+U correction for LS state. 



The natural question is rising - which spin configuration is more favorable 
- the low spin state or intermediate spin state? The energy difference (per 
Co2 ion) between these two states depended on on-site Hubbard correction 
parameter U is depicted on the Figure 5. In the range eV < U < 3 
eV it is not possible to stabilize IS state (simply, it does not exist), which 
shown as zero energy difference. For the values of U equal to 4 and 5 eV 
LS spin state has less energy, whereas for higher U value the situation is 
inverse. But neither zero point energy nor finite temperature effect was not 
taken into account, whereas, in fact, they can influence onto the ground state 
determination. 
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Table 2: Experimental lattice parameters for Na^CoC^ elementary cell. 



Compound 


Ref. 


a(A) 


c(A) 


Na . 6 iCoO 2 


[21] 


2.832 


10.843 


Na . 66 CoO2 


[22] 


2.837 


10.97 


Na .67CoO 2 


[23] 




10.938 



AE,eV 




U, eV 



Figure 5: Difference between the ground-state energy of the system in IS and 
LS states in a dependence on the parameter of on-site Hubbard interaction 
parameter U . Here, the energy is calculated per one cobalt ion in the Co2 
crystallographic position. When U is less than 4 eV, the IS state is unstable 
and only the LS state is accessible. 
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Figure 6: Band structure for LS (a) and IS (b) states about the Fermi level. 
The LS state demonstrates the metallic character, while IS shows a small 
insulating band gap (~ 30 meV). Here, the on-site Hubbard interaction pa- 
rameter U = 4 eV. 

The band structures calculated for both of the spin states are presented 
on Figure 6. IS state possess a small band gap about 30 meV at the Fermi 
level. Whereas bands structure of the LS state has the definitive metallic 
character for all the range of U. Experimentally, the metallic state of the 
cobaltate was confirmed by the measurement of the resistivity [24]. 

To calculate Fermi surface (FS) the Brillouin zone (BZ) was sampled 
with 9x9x9 k-mesh. Since band structures show no variation of energy 
in z direction, perpendicular to C0O2 plane [Z — G line on the Figure 6 
is straight), we are interested in distribution of energy in XY plane. On 
Figure 7 band structures of LS state, cut at Fermi level for two different 
values of U (U=l eV and U = 6 eV) depends on k x and k y are depicted. 
Because of the relatively rough k-spacing in Brilloiun zone, some numerical 
interpolation artefacts are arised. Nevertheless, one can note the tendency 
of increasing of volume of central Fermi pocket and decreasing of volume of 
additional small Fermi pockets, in accordance with the Luttinger theorem. 
But the six small pockets in Fermi surface remains independently to the value 
of the Hubbard U parameter. 

Moreover, while central pocket on Fermi surface always has hole char- 
acter (convex), the small side pockets change their natures from electron- 
like(concave) at U — 1 eV to hole-like (convex) at U — 6 eV, as depicted on 
Figure 7. For U=l eV, electron velocity at the central FS pocket is about 
t> F ~ 4 x 10 4 m/s, whereas for sodium concentration 1 = 1/3 «f ~ 6 x 10 4 
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m/s. The size of FS central pocket in our case is about k F ~ 0.1 A- 1 . This 
value is less compared to x ~ 1/3 case (where k-p ~ 0.7 A -1 [25]), but it 
should be noted, that the elementary cell at x = 2/3 is larger, which should 
results in decreasing of first BZ size, and thus the size of FS, deduced from 
the experiment. 

Analysis of macroscopic magnetic susceptibility of Na 2 /3Co0 2 sample in 
temperature range from 50 K to 250 K in Ref. [23] leads to effective magnetic 
moment /j, e g = 1.12/ie per Co2 site, whereas effective magnetic moment for 
LS state is only about 0.3/^b per Co2 site. One could suggest that the above- 
mentioned experimental data should correspond to IS state. But, firstly, 
these data indicate only the average magnetic moment and extra magnetic 
moments can be induced by some impurities of sample. Secondly, IS state has 
a small band gap and should not demonstrate metallic conductivity (linear 
dependence of resistance on temperature) as it observed for Nao.yCo02 [24] 
for wide temperature range. 

4 Discussion 

The existence of the strong FM correlations for the Hubbard model on the 
trigonal lattice while hole concentration corresponds to x = 2/3 has been 
shown in Ref. [26]. Such a system possess the instability with respect to 
simultaneous appearance of charge and spin density waves with periodicity 
of Kagome superlattice. The probable reason of a such type of instability is 
in van-Hoove singularity of the density of states of 3rf-electrons nearby the 
Fermi level on the Kagome lattice [27]. This explains why the extensional 
wave function of the Co 4+ hole with spin s = 1/2 tends to localize from the 
whole layer into the sublattice of Co2: it leads to decreasing of the magnetic 
energy due to the FM transition and thus, optimization of on-site Hubbard 
(Coulomb) correlations. At the same time, the inert states Co 3+ (s=0) are 
superseded to the residuary Col sites. In turn, it makes the next to the Col 
ions Nal sites attractive for the Na + ions. Then the Coulomb interaction 
between the sodium sublattice and cobalt sublattice tends to pin the whole 
charge order in the system. Such a scheme could be the probable scenario 
of stabilization in a metallic system: FM of a ground state with concurrent 
charge ordering of Co 3+ (Co 4+ ) ions on the Col (Co2) sites and arrangement 
of the Na + ions on the Nal and Na2 sites. 

Charge and spin separation for cobalt ions has been demonstrated in the 
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Figure 7: Band structure near Fermi level for LS state in XY plane for two 
cases: U=l eV (left) and U=6 eV(right). On the left figure, Fermi level 
is crossed by two bands (green and red), while on the right, only by one 
band (green). Moreover the six small pockets change its character from the 
electron-like (concave) to the hole-like (convex). 
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present work. An AFM-A type magnetic order (antiferromagnetic interplane 
and ferromagnetic inplane) has been derived as more energetically stable. We 
have studied two possible spin states - low spin state and intermediate spin 
state. Obtained IS state resemble a spin-polaron excitation state, proposed 
by Khaliullin et al. [28,29], but this state is an insulator state with a 
small band gap of about 30 meV. Our calculations have indicated that for 
small value of the Hubbard parameter U < 6 eV the ground state is a LS 
state, whereas at U = 6 eV IS state has less energy. Also, the charge/spin 
separation dependence on the Hubbard U parameter has been observed. 

The band theory predict the complicated Fermi surface (FS) with one 
large hole pocket around the Y = (0, 0, 0) point and six small pockets near the 
K = (0, -y, 0) points of the hexagonal Brillouin zone at least for x < 0.5 [9]. 
However, intensive investigations by several ARPES groups reveal absence 
of six small pockets in both Na x Co02 *yYi20 and in its parent compound 
Na^CoOs [12,24,30]. 

Independently on the value of on site Hubbard correlation energy U, the 
calculated Fermi surface has six small pockets that are not presented in the 
ARPES experiment [24]. The disagreement between ARPES spectra and ab 
initio calculated band structure points to the importance of the electronic 
correlations in these oxides. Possible cause of this unconformity is that the 
ARPES experiment shows electronic states of the sample surface, whereas, in 
certain cases, they can differ significantly from the bulk states [31]. Indeed, 
independent experiments on transport properties which exploits Shubnikov - 
de Haas effect in Na x Co02 for x = 0.71 and 0.84, exhibit the presence of the 
Fermi surface pockets [32] . Thus, for the correct understanding of contradic- 
tion between ARPES experiment and ab initio calculations of realistic crystal 
cell, an ab initio studying of sodium cobaltate surfaces should be done. 

Comparison of calculated results with the experimental data on conduc- 
tivity and Fermi surface indicates that ground state is a LS state. At the 
same time, the macroscopic magnetic susceptibility measurements indicate 
the high magnetic moment per cobalt ion and thus one can not conclude that 
one has the IS of sample. Nevertheless, it is more probable, that the extra 
magnetic moments are not appeared from Na2/3Co02 compound, but they 
can be induced by some impurities of the sample or defects of crystal lattice 
(such as sodium vacancies [10]). 
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5 Conclusion 



In the present work we have modeled electronic and magnetic properties of 
sodium cobaltate Na 2 / 3 Co0 2 within GGA+U approximation. In contrast to 
previous work [11,12,33], for the first time, a realistic, experimentally ob- 
tained elementary cell of this compound [10] was used. Ab initio calculations 
indicate that proposed sodium imprint is stable in the sense that it does not 
change after the full structure relaxation. 

The microscopic theory of FM metallic state of Na 2 /3Co0 2 with a real 
Na + arrangement in Nal and Na2 sites [10] was developed. It was shown that 
at U > 3eV in the system the complex ordering, which combines the redistri- 
bution of Co 4+ into Co2 sublattice, arises. The calculated spatial structure 
of the electronic ground state is coherent to the experimentally deduced [10]. 
It was established that the presence of both electrostatic potential induced 
by the periodic Na + arrangement and the strong Coulomb correlations in the 
3d orbitals is sufficient to appearance of such ground state. Thus, the tricky 
ground state of the cobalte could be the result of the peculiar cooperation 
between two relevant interacting subsystems: the patterning ions and the 
strongly correlated 3d electrons of the Co-layers. 
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